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An experimental study of the viscosity of a macroscopic suspension, i.e. a suspension for which 
Brownian motion can be neglected, under steady shear is presented. The suspension is prepared 
with a high packing fraction and is density-matched in a Newtonian carrier fluid. The viscosity of 
the suspension depends on the shear rate and the time of shearing. It is shown for the first time 
that a macroscopic suspension shows thixotropic viscosity, i.e. shear-thinning with a long relaxation 
time as a unique function of shear. The relaxation times show a systematic decrease with increasing 
shear rate. These relaxation times are larger when decreasing the shear rates, compared to those 
observed after increasing the shear. The time scales involved are about 10000 times larger than the 
viscous time scale r^isc = jv and about 1000 times smaller than the thermodynamic time scale 
Ttherm = P^H- The structurc of the suspension at the outer cylinder of a viscometer is monitored 
with a camera, showing the formation of a hexagonal structure. The temporal decrease of the vis- 
cosity under shear coincides with the formation of this hexagonal pattern. 

PACS number(s): 83.60.Fg, 83.60.Rs, 82.70.Kj, 47.27.Pa, 47.50. -l-d, 83.50.Ax 



I. INTRODUCTION 



The rheology of suspensions has been studied for many 
years, mainly because of its obvious importance in a wide 
range of industrial applications (see for example [H). 

A lot of work has been done especially for colloidal sus- 
pensions, i.e. suspensions where Brownian motion cannot 
be neglected i^B], §, §, §, [§, §, @, 0, @, 
p3[, 1 14 1, ]l5[|, |16[~ In these colloidal suspensions shear- 



thinning can be observed, i.e. the viscosity of the fluid 
decreases with increasing shear rate |17 , Q, [|l9|, as well 



as shear-thickening for higher shear rates. It is known, 
that colloidal suspensions also show thixotropic behavior 
pO| , [pT| , i.e. the continuous reduction in viscosity with 
time of shearing and the subsequent recovery under re- 
duced shear g, |l8|, |l|, [||, ll. The idea is proposed, 
that the change in viscosity is caused by the formation of 
structures within a colloidal suspension under shear 

, ||2l[ , . It is assumed that for colloidal suspensions 
the dynamics is dominated by Brownian diffusion [ p^ , 

MM- 

A measure for the relative importance of shear forces 
compared to Brownian motion is given by the Peclet- 
number Pe = Gir-qjr^ /{kBT), where r is the radius of 
the particles. Whereas colloidal suspensions have been 
investigated intensively up to a Peclet-number of 100, 
few experimental data exist for macroscopic suspensions, 

1. e. for suspensions of higher Peclet-number where Brow- 
nian motion can be neglected. It is assumed that for high 
Peclet numbers the shear-thinning region shifts rapidly 
to extremely low values and essentially becomes exper- 
imentally inaccessible [Q, The lack of Brownian 
motion for these large particles is believed to eliminate 
the reversibility in structural changes and should compli- 
cate experiments ||2l|] , . These difficulties explain the 
dearth of experimental data on the behavior of macro- 



scopic, i.e. non-Brownian suspensions. 

For this class interesting effects have been observed for 
non-A^ewtoman carrier fluids (see e.g. (27), H], H). 
Suspensions in Newtonian media appear to have been 
almost neglected. Some works showed the difficulties 
arising when trying to investigate the behavior of con- 
centrated suspensions. It was demonstrated that shear- 
ing the macroscopic suspensions in a concentric cylinder- 
viscometer of the Searle-type leads to shear migration ef- 
fects, i.e. the particles migrate in nonhomogeneous shear 
flows from high shear rates to low shear rates [Q, 
p2[ p3[ . For the geometry used in those experiments it 
meant that the particles even migrated out of the sheared 
volume and the particle concentration decreased. Since 
for high packing fractions the viscosity depends strongly 
on the concentration, this reduction in concentration was 
measured. So from these measurements no estimate of 
the real viscosity of the suspension could be made. 

In another experiment a macroscopic suspension 
was exposed to oscillatory shear and the decrease of the 
viscosity with time was measured. In this experiment 
there was a hint for an increase in the organization of 
the particles with shearing time. Nevertheless the sus- 
pension was not density matched, so it is not proved that 
the decrease in viscosity is not, at least partially, due to 
sedimentation. 

The present paper presents experimental results on the 
behavior of macroscopic suspensions under steady shear 
at 7 = 50/s of a Peclet-number Pe « 2.1 • 10^ when us- 
ing the viscosity of the suspension or Pe « 8.3 • 10^ when 
using the viscosity of the carrier fluid, i.e. the Peclet num- 
ber is seven orders of magnitude larger than the Peclet 
numbers which are covered by the experiments with col- 
loidal suspensions. We flnd evidence for a thixotropic 
behavior of the suspension with extremely large relax- 
ation times. The relaxation times are about 10000 times 
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larger than the viscous time scale t^sc = jv and about 
1000 times smaller than expected from the Peclet number 
Ttherm — PP'/j- This makcs it unlikely that a hydrody- 
namic or thermodynamic interaction is responsible for 
the effect. One can thus speculate that the time scale is 
dominated by ordering effects within the densely packed 
granular network provided by the beads |^^. The de- 
crease of the viscosity under shear coincides with an in- 
crease in order of the particles. 

The paper is organized as follows: In the following sec- 
tion the experimental setup is described, in Section III 
the experimental results are discussed. A final section 
contains the conclusions. 
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FIG. 1. Size distribution of the used particles, measured 
from conductivity measurements. 



II. EXPERIMENTAL SETUP 



The suspension investigated consisted of nearly 
monodisperse glass beads (Wiirth Ballotini MGL) dis- 
persed in an aqueous solution of zinc iodide Znl2. The 
glass beads are spherical and consist mainly of 5*102 and 
have a material density of p = 2.45 g/cm"^. The parti- 
cles are sieved in order to obtain a well defined particle 
size distribution which is then measured with a Multi- 
sizer 2 (Coulter). The particle size distribution (PSD) is 
shown in Fig. |[ This measurement of the PSD is based 
on the principle of monitoring the electrical current be- 
tween two electrodes immersed in a conductive liquid on 
either side of a small aperture, through which the sus- 
pension of particles is forced to flow. As a particle passes 
through the aperture, it changes the impedance between 
the electrodes and produces an electrical pulse of short 
duration with an amplitude essentially proportional to 
the particle volume. The arithmetic mean of the particle 
size is 61.13 fim and the standard deviation is 18.71 /im. 
An alternative measurement of the particle size distri- 
bution was based on laser light diffraction experiments. 
These measurements confirm the results in Fig. |^. Before 
the glass beads were suspended in the solution they were 
washed 8 times in totally demineralized water, in order 
to purify them from any residues and dried afterwards. 



To stabilize the suspension, i.e. to avoid sedimenta- 
tion during the time-consuming rheological experiments, 
the density of the continuous phase has been matched 
to that of the glass beads. Hence the particles were sus- 
pended in almost saturated Zn/2-solution, which had the 
same density p = 2.45 g/cm^ at T = 20.0°C as the glass 
beads. (The maximal density of the solution is about 2.6 
g/cm'^.) The density-matched, i.e. stabilized, suspension 
was prepared with the packing fraction cj) = 50.49%. It 
was checked at T = 20.0°C for several hours that the sus- 
pension stayed homogeneous and the glass beads did not 
sediment to the bottom or rise to the top of the vessel. 
The carrier fluid, the pure Z?i/2-solution, is a Newtonian 
solution with a dynamic viscosity of 12 mPas at 20.0°C. 

The rheological measurements were performed with a 
rheometer SR-5000 (Rheometric Scientific, Inc.) with a 
rotating inner cylinder in a fixed outer cylinder. The in- 
ner cylinder has a cone at the lower end. Due to the 
relatively large particles a concentric cylinder geometry 
with a ratio 5 — Ri/Ro — 0.922 of the inner cylin- 
der radius Ri — 14.75 mm to the outer cylinder radius 
Ro = 16.0 mm and a gap width of a=1.25 mm has been 
used. The length of the inner cylinder is 44.3 mm. The 
shear stress in the small gap is constant within 15% for 
a Newtonian fiuid. The thixotropic behavior was inves- 
tigated in the controlled shear rate mode by a stepwise 
increase and subsequently stepwise decrease of the shear 
rate 7 and registration of the shear stress a and the ap- 
parent viscosity rjapp — cr/j. The temperature during the 
whole measurements was stabilized to 20.0 °C. 

The shear rates are below 7 ~ 500/s, which corre- 
sponds to a Taylor number |||] Ta = n^Ra^v-^ « 11.3, 
where ft is the angular velocity and R = {Ri + Ro)/2. 
Hence the experiments are performed far away from the 
onset of Taylor-vortices at Ta=1712. 

In order to allow for an optical investigation of the 
structure-formation at the outer cylinder of the suspen- 
sion, a shearing-apparatus was built which had the same 
gap width a=1.25 mm and almost the same ratio d— 0.916 
as the rheometer. The outer cylinder had a radius of 15.0 
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mm and was made out of glass to allow the view into the 
suspension layer at the outer cylinder. The inner cylinder 
did not have a conic but a flat lower end. Images were 
taken during shearing. In addition the torque of the mo- 
tor was monitored, which is a measure for the viscosity 
of the suspension. 



III. EXPERIMENTAL RESULTS 
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FIG. 2. Thixotropic behavior of a suspension of packing 
fraction of 50.49%. (a) Driven shear rate ramp, (b) stress, 
(c) dynamic viscosity of the suspension 



The following sequence of shear rate steps was used: 
7 = 0/s, 5/s, 50/s, 5/s, 0/s. The resulting shear stresses 
and the apparent viscosities rjapp were measured as a 
function of time. Fig. |^(a) shows the shear rate steps 
together with the system response, i.e. (b) the stress and 
(c) the dynamic viscosity rj. When 7 is set from to 5/s 
an exponential decay of the viscosity to the final viscos- 
ity is observed. Then the shear rate is set to 50/s and 
another exponential decay to another apparent viscosity 
is observed. This final value of the viscosity is smaller 
than the final value at 7 =5/s. After setting 7 back to 
5/s the former final viscosity at this shear rate is gradu- 
ally reached again. This is the behavior of a thixotropic 
material. 



From this measurement the relaxation times r are ob- 
tained from a fit to 



■nit) = a- exp{-t/T) +7]fi„ai (3.1) 



Fig. II and Fig. ^ show the fits for 7 =50/s (from 
7 =5/s) and for 7 =5/s (from 7 =50/s). 
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FIG. 3. Exponential fit to tlie decay of the dynamic vis- 
cosity with time for the step from 7 =5/s to 7 =50/s. For 
the reason of clearness only every 10th data point is shown. 
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FIG. 4. Exponential fit to the growth of the dynamic vis- 
cosity with time for the step from 7 =50/s to 7 =5/s. For 
the reason of clearness only every 10th data point is shown. 

For the first shear rate 7 =5/s the relaxation time is 
T=1276 s. For the shear rate 7 =50/s the relaxation time 
is T=968 s and for the gradual rise in viscosity for 7=5/s 
the relaxation time r=12182 s, i.e. it is one order of mag- 
nitude larger than for the decay in viscosity at the same 
shear rate. 

The ramp for these shear rates was not applied in a 
repeated sequence with the same suspension, but was ap- 
plied four times for freshly prepared suspensions and the 
thixotropic behavior of the suspensions was reproducible. 
For the following shear rates the ramp was applied once 
with alternating jumps from higher to lower and lower to 
higher shear rates. 

This evaluation of the relaxation times was also done 
for other measurements for both possible cases: coming 
to the appropriate shear rate from a lower shear rate (step 
from 7=0 to 5/s, from 5/s to 50/s and from 10/s to 100/s) 
or coming from a higher shear rate (step from 7=100/s to 
50/s and 500/s to 100/s). The duration of shearing was 



3 hours when coming from a lower shear rate and 6 hours 
when coming from a higher shear rate, respectively. The 
relaxation times show a systematic decrease for increas- 
ing shear rates. For decreasing steps from higher shear 
rates to lower shear rates the relaxation times are larger 
than for increasing steps. All results are shown in Fig. |^. 
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FIG. 5. Dependence of the relaxation time r on the shear 
rate. Filled symbols represent steps from higher shear rates 
to lower shear rates, open symbols represent steps from lower 
shear rates to higher shear rates. 

Using the optical setup, the torque of the motor was 
monitored. Fig. ^ shows the decrease of the torque aver- 
aged over 100 s. An exponential fit according to equ. 3.1 
gives the relaxation time of 1089 s, which is close to the 
one measured with the viscometer. Although the sig- 
nal to noise ratio is worse than with the viscometer this 
allows the interpretation that the same mechanism is re- 
sponsible for the torque reduction in both apparatuses. 
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FIG. 6. Decrease of the torque of the motor, which is a 
measure of the viscosity of the suspension under shear. The 
torque is averaged over 100 s. Every 3rd data point is shown. 

Images with a Long Range Microscope (Questar 
QMIOO) which is connected to a CCD-camera using an 
exposure time of about 1 ms are taken from the suspen- 
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sion at the outer cylinder every 10 s during shearing at 
7 =50/s. The starting time of the shearing is defined 
by t = 0. The images show a cantle of 1.53 x 1.53 mm 
and are enhanced in contrast. Fig. 0(a) shows images 
of the suspension at the outer cyhnder at t = and at 
t = 62 min. There is almost no difference observable 
in real space but the 2D-Fourier spectra of the images 
(Fig. 0(b)) show a hexagonal structure after 62 min of 
shearing. 
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FIG. 7. (a) Images of the suspension at the outer cylinder at t = and t — 62 min. (b) 2D-Fourier spectra of the 
images. At t = no structure is exhibited, at t = 62 min a hexagonal structure is dominant. The shear is applied in the 
a;-direction and the axis of the rotating cylinder is parallel to y. 
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The increase of this hexagonal structure is quanti- 
fied by calculating a measure for the contrast a as the 
root-mean-square normalized by the mean intensity of 

the spectral intensities for wave numbers k — ^ k"^ + fc^ 

within a range of 94/mm < k < 106/mm, i.e. 

i V " 

where i is the intensity in the spectrum and n is the 
number of data points. This contrast is averaged over 20 
images in order to reduce the noise. The unstructured 
spectra have a lower contrast a than the structured spec- 
tra, i.e. a is a measure for the structure of the suspension 
at the outer wall of the cylinder. Fig. ^ shows the tem- 
poral increase of that number. 




IV. CONCLUSION 



An experimental study of the viscosity of a macro- 
scopic suspension of high packing fraction which is 
density-matched in a Newtonian carrier fluid under steady 
shear is presented. It is demonstrated that macroscopic 
suspensions can be considered as thixotropic materials. 

Furthermore relaxation times at the order of hours 
were observed. They show a systematic decrease for in- 
creasing shear rates. For decreasing steps from higher 
shear rates to lower shear rates the relaxation times are 
larger than for increasing steps. 

During the decrease of the viscosity under shear a 
structure-formation was observed. The Fourier spectra 
show that the system forms a hexagonal structure within 
time. Both effects, i.e. the decrease of viscosity and the 
increase of patterns, coincide. 

We cannot provide an explanation for the physical ori- 
gin of the thixotropic behavior. We note that the relax- 
ation times are about 10000 times larger than the viscous 
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time scale t^sc = d^jv and about 1000 times smaller 
than the thermodynamic time scale Ttherm = P^/j- This 
makes it unlikely that hydrodynamic or thermodynamic 
interactions are responsible for the effect. 

One can thus speculate that the time scale is domi- 
nated by ordering effects within the densely packed gran- 
ular network provided by the beads p5|| . In particular, 
we believe that the friction between the wetted surfaces 
of spheres is primarily responsible for the measured ef- 
fective viscosity. In a densely packed granulate material 
that friction will be dominated by the fact that individ- 
ual particles are subject to frustration if they feel more 
than one particle contact, which tends to increase the ef- 
fective viscosity. Ordering under shear might mean that 
the particles reduce the mutual friction by reducing the 
number of frustrated contacts which in turn might cause 
a lower viscosity. This idea is related to the concept 
of dilatancy, which is fundamental in the description of 
fluid properties of granular matter. Additional investi- 
gations using three-dimensional imaging techniques and 
theoretical calculations are clearly necessary to substan- 
tiate these speculations. 
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